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ABSTRACT 

In recent years, a considerable research effort has been devoted to the study of heat transfer induced by natural 
convection within a cavity. 

The interest in these phenomena of natural convection is due to the numerous potential applications in 
engineering. These applications include the extraction of geothermal energy, dispersion of pollutants in aquifers, safety 
problems in the core of nuclear reactors, and thermal building. In recent years, natural ventilation, similar to 
the phenomenon of the natural convection of living quarters, is approached numerically using primarily the CFD code. It is 
in this context that our study is located and wants an additional contribution and which is centered around the modeling 
and the simulation by finite volumes of a system of natural ventilation of a room in a region arid by a composite wall Used 
in the ventilation system or the passive heating of buildings (analysis of the natural convection phenomena). 

The results are presented in the form of temperature fields and current lines for Rayleigh number 10 5 , air flow 
mass and local and average Nusselt numbers. 
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MODELING OF NATUREL CONVECTION 

Convection in two-dimensional cavities (rectangular or square), differentially heated by their vertical faces 
(Window Problem or Double Glazing) has been studied extensively over the past twenty years. 

Modeling Assumptions 

Our modeling is based on the following simplifying assumptions: 

• The fluid is Newtonian; 

• The fluid is incompressible; 

• The regime is laminar; 

• Heat transfer is two-dimensional; 

• The thermo-physical properties of the fluid are constant; 

• The dissipation of energy in viscous form is neglected. 

• The fluid satisfying the Bossiness approximation. 

p(T) = p„{l-P T (T-T„)\ 
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The decomposition of conservation equations into Cartesian coordinates gives the following system of differential 
equations: 

Conservation of the Mass 


du dv 
— + — = 0 

dx dy 


Conservation of the Quantity of Movement 
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Conservation of Energy 
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In these equations, u and v are the components of the velocity, P the pressure, T the temperature, to the time, g the 
acceleration component of gravity, the thermal diffusivity and the density 

By choosing characteristic scales for each of the variables of the problem, these equations can be dimensioned, 
which gives dimensionless numbers characteristic of natural convection (. Ra , PR...). 

The system of equations governing the flow of the fluid is then written in the non-dimensional form: 


dir av* 
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= 0 
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All the equations that govern our problem can be written in the following conservative form: 
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A 

dt 


f- 


V.V 
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= v 




+ ^ 


/ \ 

Where (j) the variables of the problem represent \U , V , T J , and S ^ Are the diffusion coefficient and 
the source term, respectively. 

CHOICE OF THE MESH 


Le maillage utilize est. un maillage uniform de 61*61 nudes, suivant les directions x et y. 



— = 0 , u - v — 0 

ay 



Figure 1: Geometry and Boundary Conditions 


CONDITIONS TO THE LIMITS IN THE OPENINGS OF THE ENTRY AND THE EXIT 

In the pressure inlet boundary state, Fluent used to define the fluid pressure at the inlet flow all other scalar 
properties of the flow. 

When the inlet pressure is known, but the flow rate or speed is not known. Fluent uses a technique to define a free 
boundary condition; this situation can arise in many practical situations. 
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In the incompressible flow, the total inlet pressure P 0 and the static pressure P s , are related to the rate of 


admission by the Bernoulli equation 


Po=Ps+^PV : 


For incompressible flows, the boundary pressure is calculated as an average between the external pressure and the 
internal pressure: 

P, = 0.5 (P c +P.) 

METHOD OF RESOLUTION 


Numerical resolution of partial differential equations has developed considerably since 1940 with the 
development of computers that are increasingly fast and have a larger memory size. 


A large majority of commercial computation codes solves the conservation equations in primitive variables and 
use the finite volume or control volume method. 


In this method, the domain under consideration is divided into finite volumes; a single volume surrounds a mesh 
node. The conservation equations are integrated on each finite volume. The model thus discretized, will be solved under 
the CFD (Fluent) code. 

The simulation was carried out in transitional regimes in order to study the natural ventilation in a cavity with a 
composite wall. 

The parameters of the simulation are as follows: PR = 0.71 (air), Ra = 10 5 

The initial temperature in the cavity TO = 0 (see FIG. 2). 
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Figure 2: Initial Temperature At (T = Os) 
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Temperature Profiles 



Figure 4: Profiles of Temperature T in the Horizontal 
Median Plane for Different Numbers of Ra 



Figure 5: Profiles of Temperature T in the Vertical 
Median Plane for Different Numbers of Ra 
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Variation of air flow at Entrance 

For the basic situation, we present the mass of airflow as a function of the Rayleigh number in Figure 6; we see 
that the mass of air entered with the increase of the Rayleigh number. 

These values show that the air suction increases with the increase of Rayleigh number, so the mass of air flow 
tends to a maximum value. 



Figure 6: Variation of the Mass of Air Flow at the 
Entrance for Different Numbers of Ra 


Variation in Local Nusselt Number 

La figure 7 shows the evolution of the heat transfer (local Nusselt number) as a function of y * at the hot wall of 
the wall for different Rayleigh numbers Ra. 



Figure 7: Local Nusselt Number for Different Number of Ra 

Variation in Number of Medium Nusselt 

Figure 8 shows the average Nusselt number as a function of Rayleigh number, these values show that the average 
Nusselt number increases with the increasing Rayleigh number, and tends to a maximum value. Heat transfer by natural 
convection becomes more dominant, so the average Nusselt number is small forward 
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Low Rayleigh number (10 3 - 10 4 - 10 5 ), for higher Rayleigh numbers 10 6 the average Nusselt number becomes 

higher. 

These observations are expected the fireplace system to become perform better when the Rayleigh number 
increases, which indicated heat transfer and ventilation 



Figure 8: Variation of the Average Nusselt 
Number for Different Number of Ra 


CONCLUSIONS 

Based on this model, we have reached the following conclusions: 

The effect of the Rayleigh number that characterized the intensity of natural convection was highlighted. Thermal 
transfer increases with the Rayleigh number. 

In the study that was presented, the natural ventilation system seems like a very interesting solution for: 

• Thermal comfort and olfactory comfort. 

• Energy saving (passive solar energy). 

• To eliminate the pollutants (odors, water vapor, gas,). 

• To remove heat. 
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